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We construct systematic expansions around four and two spatial dimensions for a Fermi gas near 
the unitarity limit. Near four spatial dimensions such a Fermi gas can be understood as a weakly 
interacting system of fermionic and bosonic degrees of freedom. To the leading and next-to-leading 
orders in the expansion over e = 4 — d, with d being the dimensionality of space, we calculate the 
thermodynamic quantities and the fermion quasiparticle spectrum, both at unitarity and around the 
unitarity point. Then the phase structure of the polarized Fermi gas in the unitary regime is studied 
in the e expansion. At unitarity the unpolarized superfluid state and a fully polarized normal state 
are separated by a first order phase transition. However, on the BEC side of the unitarity point, in a 
certain range of the two-body binding energy, these two phases are separated in the phase diagram 
by more exotic phases, including gapless superfluid phases with one or two Fermi surfaces and a 
superfluid phase with a spatially varying condensate. We also show that the unitary Fermi gas near 
two spatial dimensions is weakly interacting, and calculate the thermodynamic quantities and the 
fermion quasiparticle spectrum in the expansion over e — d — 2. 

PACS numbers: 03.75.Ss, 05.30.Fk 
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I. INTRODUCTION 

Two-component Fermi gas with zero-range interaction 
at infinite scattering length [1, 2], frequently referred to 
as the unitary Fermi gas, has attracted intense attention 
across many subfields of physics. Experimentally, the 
system can be realized in atomic traps using the Fesh- 
bach resonance and has been extensively studied [3-9]. 
Since the fermion density is the only dimensionful scale 
of the unitary Fermi gas, its properties are universal, 
i.e., independent of details of the interparticle interac- 
tion. The unitary Fermi gas is an idealization of dilute 
nuclear matter and may be relevant to the physics of 
neutron stars [10]. It has been also suggested that its 
understanding may be important for high-Tc supercon- 
ductivity [11]. 

The austere simplicity of the unitary Fermi gas im- 
plies great difficulties for theoretical treatment, because 
there seems to be no parameter for a perturbation the- 
ory. The usual Green's function techniques for the many- 
body problem are completely unreliable here since the 
expansion parameter akp becomes infinite in the uni- 
tarity limit. Considerable progress has been made by 
Monte Carlo simulations [12-15], however these simula- 
tions have many limitations, which become especially ev- 
ident for systems with a population imbalance between 
two fermion components (finite polarization). 

Recently, we have proposed a new approach for the 
unitary Fermi gas based on the systematic expansion 
in terms of the dimensionality of space [16], utilizing 
the specialty of four spatial dimensions in the unitar- 
ity limit [17]. In this approach, one would extend the 
problem to arbitrary spatial dimension d and consider d 
is close to, but below, four. Then one performs all cal- 
culations treating e = 4 — d as a small parameter of the 
perturbative expansion. Results for the physical case of 



three spatial dimensions are obtained by extrapolating 
the series expansions to e = 1. 

We used this e expansion around four spatial dimen- 
sions to calculate thermodynamic quantities and fermion 
quasiparticle spectrum in the unitarity limit and found 
results which are quite consistent with those obtained 
by Monte Carlo simulations and experiments [16]. Very 
recently, this expansion has been successfully applied to 
atom-dimer and dimer-dimer scatterings in vacuum [18]. 
Thus there are compelling reasons to hope that the limit 
d — > 4 is not only theoretically interesting but also prac- 
tically useful, despite the fact that the expansion param- 
eter e is one at d = 3. One may hope that the expansion 
over 4 — d will be as fruitful for the Fermi gas at unitar- 
ity as it has been in the theory of the second order phase 
transition [19]. 

This paper has four main purposes. We will give a 
detailed account of the e expansion for the unitary Fermi 
gas studied in our recent Letter [16], extend the results 
to the Fermi gas with large but finite scattering length, 
study a phase diagram of a polarized Fermi gas near the 
unitarity limit, and develop a complementary expansion 
for the unitary Fermi gas around two spatial dimensions. 

The structure of the paper is as follows. We start with 
the study of two-body scattering in vacuum in the unitar- 
ity limit for arbitrary spatial dimension 2 < d < 4. This 
study will clarify why systematic expansions around four 
and two spatial dimensions are possible for the unitary 
Fermi gas (Sec. II). The e expansion is developed for an 
unpolarized Fermi gas both at unitarity and at large, but 
finite, scattering length in Sec. III. Then we apply the e 
expansion to the unitary Fermi gas with unequal densi- 
ties for the two fermion components. Such a system will 
be called "polarized" Fermi gas and has been recently re- 
alized in experiments [20-23] . The phase structure of the 
polarized Fermi gas in the unitary regime is investigated 
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based on the e expansion in Sec. IV. We also show in 
Sec. V that there exists a systematic expansion for the 
Fermi gas in the unitarity hmit around two spatial di- 
mensions. In Sec. VI, we make an exploratory discussion 
to connect the two expansions around four and two spa- 
tial dimensions. The summary and concluding remarks 
are given in Sec. VII. 



II. TWO-BODY SCATTERING IN VACUUM 
A. Qualitative discussion 

Nussinov and Nussinov [17] were the first two to rec- 
ognize the special role of two and four spatial dimensions 
in the unitarity limit. Their arguments are very simple 
and we review them here for completeness. 

In low dimensions d < 2, any attractive potential pos- 
sesses at least one bound state. Therefore, the threshold 
of the appearance of the first two-body bound state cor- 
responds to the zero coupling. It follows that the Fermi 
gas in the unitary limit corresponds to a noninteract- 
ing Fermi gas at d ^ 2 and the energy per particle ap- 
proaches that of the free Fermi gas in this limit. If one 
defines a parameter ^ as the ratio of the energy density 
of the Fermi gas at miitarity and that of the free Fermi 
gas with the same density. 
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(1) 



then according to the above argument ^ — > 1 as d ^ 2 
from above. (The singular character of d = 2 was also 
recognized in the earlier work [24].) 

Now we consider dimensions close to four. At infi- 
nite scattering length, the wave function of two fermions 
with opposite spin behaves like R{r) = l/r'^~'^ at small r, 
where r is the separation between two fermions. There- 
fore, the normalization integral of the wave function takes 
the form 



j drR{rf = ^ 



1 



,d-3 • 



(2) 



which has a singularity at r ^ in high dimensions d > 
4. For these dimensions the two-body wave function has 
infinite probability weight at the origin, and the ferniion 
pair looks like a pointlike boson. From this observation, 
Nussinov and Nussinov concluded that the unitary Fermi 
gas at (i ^ 4 becomes a noninteracting Bose gas. In 
particular, the energy per particle at fixed Fermi energy 
goes to zero as d ^ 4, because all bosons condense into 
the zero energy state and their binding energy is also 
zero. According to this argument, ^ ^ as d — > 4 from 
below. 

This qualitative observation by Nussinov and Nussi- 
nov, however, leaves many questions unanswered: 

1. How do the bosons interact with each other? 



2. How fast does ^ approach as d — > 4? 

3. Can one develop a theory based on an expansion 

over e = 4 — d? 

In order to develop a framework that can be used to an- 
swer these questions, we need to review the two-body 
scattering in a general number of spatial dimension d, 
with special attention to d near four and two, from the 
perspective of the diagrammatic approach. This consid- 
eration will provide a foundation for the construction of 
systematic expansions for the many-body problem of the 
unitary Fermi gas. 



B. Around four spatial dimensions 

One can describe the system of spin-i fermions with 
short-range interaction by the following Lagrangian den- 
sity (here and below h=l): 



C 



2m 



V'a + co^|f/'|V'iV'T- (3) 



Here m is the fermion mass and cq is the bare attractive 
coupling between two fermions. The last term in Eq. (3) 
corresponds to the interaction potential 



V{r) = -co6{r). 



(4) 



It is well known for d = 3 that in order to have a fi- 
nite scattering length, the potential in Eq. (4) has to be 
regularized: One chooses a particular shape of the poten- 
tial (for example, the square- well shape), and then takes 
the width of the potential to be zero simultaneously with 
adjusting the height so that the zero-energy wave func- 
tion approaches a finite limit. We shall show that it is 
true for d between two and four considering the two-body 
scattering. 

The T-matrix of the two-body scattering is given by 
the infinite summation of bubble diagrams. As a result, 
its inverse is 

1 1 , • /" dkodk 1 



=--/ 

Co Jk 



(5) 



where Sp = /2m is the kinetic energy of nonrelativistic 



particles and we introduced the shorthand notation 

r _ dk 



(6) 



The integration over k in Eq. (5) is ultraviolet divergent 
as /dfe//c^ if d > 2 and needs to be regularized. 
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The usual way to do this regularization is to assume 
that there is an upper limit A in the momentum integral 
of Eq. (5), and to adjust the A dependence of cq so that 
the physics does not depend on A. If we add and subtract 
the value of the integral at po = P = 0, we have 



1 



CO (A) 

A 



1 

fe 2efc 



(7) 



1 



The second integral is now Jdk/k'^ at large k and con- 
verges for d < 4. Thus for d < 4 one can set the cutoff to 
infinity in the second integral. The first integral is still 
ultraviolet divergent but is independent of po and p, and 
it can be absorbed into the definition of the renormalized 
coupling c"" as follows: 





FIG. 1: Two-fermion scattering in vacuum in the unitarity 
limit. The T- matrix near four spatial dimensions is expressed 
by the propagation of the boson with the small effective cou- 
pling g, while it reduces to a contact interaction with the 
small effective coupling near two spatial dimensions. 
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The zero range limit corresponds to sending A — > cx3 but 
keeping the renormalized coupling Cq°" finite by an ap- 
propriate choice of the bare coupling co(A). In this limit 
all physical quantities should depend only on Cg''" but 
not on Co and A separately. From the above discussion 
it is clear that this limit can be taken only for c? < 4; for 
d > 4 the divergences in the fermion loop can no longer 
be compensated by the renormalization of one coupling 
constant cq. There is no universal zero range and large 
scattering length limit for d> 4. 

In this paper we will not use the momentum-cutoff reg- 
ularization described above. Rather we shall use the di- 
mensional regularization, in which the integrals are eval- 
uated in dimensions where they are ultraviolet finite and 
then the results are analytically continued to the dimen- 
sion d. For example, the integral 



(9) 



converges only for d 
analytical formula 



2m < 2n, but we will use the 



r(2 



2r(i)r 



27r 



d/2 



(10) 



for all values of d. This regularization scheme is partic- 
ularly convenient because the divergent integral in the 
formula defining the renormalized coupling Cq"'^ is zero. 



1 

k 2efc 



in dimensional regularization, 



(11) 



and hence the renormalized coupling is just the bare cou- 
pling, Co°" = Cq. In other words, in the dimensional regu- 
larization the subtraction of the leading divergence of the 



fermion bubble diagram is implemented automatically. 
As long as physical quantities are insensitive to the ul- 
traviolet cutoff, which is the case for d < 4, the results 
of the dimensional and momentum-cutoff regularizations 
coincide. 

In the dimensional regularization, the integral in 
Eq. (5) vanishes for po ~ p — (scatterings at threshold). 
Therefore the limit of infinite scattering length where 
T(0, 0) = oo corresponds to cq — oo. In this limit the 
integration over k in Eq. (5) can be evaluated explicitly. 



Tipo^p)- 



-r 1 



i6 



©'(--? 

(12) 

This expression has a pole at d = 4 originating from the 
ultraviolet divergence of the k integration. Thus the T- 
matrix is small when d is close to four and vanishes at 
d = 4, indicating that such a limit can be thought of 
as a noninteracting limit. Substituting d = 4 — e and 
expanding in terms of e, the T-matrix near four spatial 
dimensions becomes 



^Tipo,p) 



m? pq 



+ 0(6^). 



(13) 



From this equation it is clear what is special about d 
close to four: As the number of dimensions approaches 
four from below, the two-fermion scattering amplitude 
has the form of the scattering process that occurs through 
the propagation of an intermediate boson whose mass is 
2m, as depicted in Fig. 1. 

The propagator of the intermediate boson is 



D 



{Po,p) = (po 



y 



iS 



(14) 



and the effective coupling of the intermediate boson to 
the fermion pair is 



(15) 
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Then the T-matrix, to leading order in e, can be written 

as 



iT{pQ,p) ~ {igfiD{po,p). 



(16) 



It is natural to interpret the intermediate boson as a 
bound state of the two fermions at threshold, which will 
be referred to simply as a boson. 

An important point is that the effective coupling of 
the two fermions into the boson q ^ small near 

four dimensions. This fact indicates the possibility to 
construct a perturbative expansion for the unitary Fermi 
gas near four spatial dimensions in terms of the small 
parameter e. 



C. Around two spatial dimensions 

Similarly, a perturbative expansion around two spatial 

dimensions is possible. This is because the inverse of 
the T-matrix in Eq. (12) has another pole at d = 2 and 
hence the T-matrix vanishes in this limit, which indicates 
that d = 2 from above is again a noninteracting limit. 
Substituting d = 2 + e and expanding in terms of e, the 
T-matrix near two spatial dimensions becomes 



iT{po,p)=i—e + 0{e^). 
m 



(17) 



If we define the effective coupling at 2 -|- e dimensions as 



2 27r 

r = — e, 

m 



(18) 



then the T-matrix to the leading order in e can be written 
as 



?:t(po:P) ^ ig^ 



(19) 



We see that the T-matrix near two spatial dimensions 
corresponds to a contact interaction with the small ef- 
fective coupling ^ e, as depicted in Fig. 1. In this 
case, the boson propagator D{p) in Eq. (16) is just a 
constant —1. We defer our discussion of the expansion 
over e = rf — 2 to Sec. V and concentrate on the expansion 
over e = 4 — d. 



D. Binding energy of two-body state 

We shall be interested not only in the physics right 
at the unitarity point, but also in the vicinity of it. In 
other words. I/cq can be nonzero in the dimensional reg- 
ularization. The case of co < corresponds to the BEC 
side of the unitarity point, and cq > corresponds to the 
BCS side. In order to facilitate an extrapolation to three 
spatial dimensions, one would like to identify cq with the 
scattering length. However, this cannot be done directly 
since the dimensionality of cq changes with the number 
of dimensions. 



One way to circumvent the problem is to parametrize 
the deviation from the unitarity by the binding energy 
of the two-body state. Clearly, this method works only 
on the BEC side of the unitarity point. Fortunately, as 
we will find in Sec. IV, all interesting exotic phases of 
the polarized Fermi gas appear on the BEC side of the 
unitarity limit. 

The binding energy eb, defined to be positive Eb > 0, 
is obtained from the location of the pole of the T-matrix 
at zero external momentum: T(— £b, 0)~^ = 0. From 
Eqs. (5) and (12), we see the coupling cq < is related 
with £b via the following equation: 



= r 1 



Co 



(20) 



Near four spatial dimensions, the relationship between Sh 
and Co to the leading order in e = 4 — is given by 



Co ~ 2e V27r/ 



(21) 



In three spatial dimensions, the relationship between 
and the scattering length a = — mco/47r becomes = 

At finite density in general d dimensions, we shall use 
Eb/^F as the dimensionless parameter characterizing the 
deviation from the unitarity, where ep is the Fermi en- 
ergy. If we recall that the Fermi energy is related to the 
Fermi momentum fcp through ep = fcp/(2m), one finds 
the following relationship holds at = 3: 



£b 



(22) 



We will use these relations to compare results in the ex- 
pansion over e = 4 — d with the physics in three spatial 
dimensions. 



III. FERMI GAS NEAR UNITARITY 

A. Lagrangian and Feynman rules 

Employing the idea of the previous section, we now 
construct the e expansion for the Fermi gas near the uni- 
tarity limit. We start with the Lagrangian density in 
Eq. (3) and introduce two chemical potentials /if and ni 
for the two spin components as follows: 

C= xPl(idt + ^+ Ma) + co^|^|^iVT- (23) 

After making a Hubbard-Stratonovich transformation, 
we rewrite the Lagrangian density as 



- —(jfcj) + -I- 

Co 



(24) 
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where ^' = (■0|,V'|)'^ is a two-component Nambu- 
Gor'kov field, and (Ji^2.3 and a± — \{cri ± i(j2) are the 
Pauh matrices. We define the average chemical potential 
as /i = + /ix)/2 and the chemical potential difference 
a.s H = {fi^ - fii)/2. 

The ground state at the finite density system (at least 
when _ff = 0) is a superfluid state where 4> condenses: 
((/)) = 00 with 00 being chosen to be real. With that in 
mind we expand </> around its vacuum expectation value 
00 as 
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7r2e)l/2 / 77100 ^'^^ 



2n 



(25) 



Here we introduced the effective coupling g ~ e^/^ in 
Eq. (15). The extra factor (7)7,0o/27r)'^^ was chosen so 
that if has the correct dimension of a nonrelativistic 
field [39]. 

We now split the Lagrangian into three parts, C = 
Co + Ci + C2, where 



2777 

4777 



CT_0O ^ 



</J 

Co 



(26) 



^2 ^ 



2^ 



9_ 

Co 



Co 



■(^ • 

Co 



4771 



(27) 



(28) 



The Lagrangian density in Eq. (24) does not have the ki- 
netic term for the boson field p. We add such a term to 
£0 and subtract it in £2- Analogously, we add a chemi- 
cal potential term for in Lx and subtract it in £2- At 
first sight, the split (26)-(28) may appear counterintu- 
itive, but it will greatly simplify the organization of the 
expansion in powers of e. 

We shall also see that the condensate 0o coincides, to 
the leading order in e, with the energy gap in the fermion 
quasiparticle spectrum. From Eq. (21), —g^/cQ ~ Cb 
gives the binding energy of the boson to the leading order 
in e when co is negative. Throughout this paper, we 
consider the vicinity of the unitarity point where £b ^ 
£00- 

The part £0 is the Lagrangian density of noninteract- 
ing fermion quasiparticles and bosons with the mass 277?.. 
The propagators of fermion and boson are generated by 
£0. The fermion propagator G is a 2 x 2 matrix. 



G(po,p) 



£;2 



i5 



Po 



-<Po 
H -e. 



(29) 



iG 



iD 



X ■- 



- ^ - 



-Z50o/co 



-i30o/co 



FIG. 2: Feynman rules from the Lagrangian density in 
Eqs. (26)-(28). Tlie two vertices in the last column come 
from £2, while the rest come from £1. Solid (dotted) lines 
represent the fermion (boson) propagator iG (iD). 



where Ep = 



; + 0Q is the excitation energy of the 

fermion quasiparticle. The boson propagator D is given 
by the same form as in Eq. (14), 



Dipo,p) = (po - y + «<5 



(30) 



The part £1 generates the vertices depicted in the first 
two columns of Fig. 2. The first two terms in £1 describe 
the fermion-boson interactions, whose coupling is propor- 
tional to g ~ e^/^ and small in the limit e ^ 0. This cou- 
pling originates from the two-body scattering in vacuum 
in the unitarity limit which was studied in the previous 
section (see Fig. 1). The third and fourth terms are in- 
sertions of chemical potentials, /i and /xb = 2/i — /cq, 
to the fermion and boson propagators, respectively. We 
treat these two insertions as small perturbations for the 
following reasons. First, we shall see that /i is small com- 
pared to the scale set by 0o, i.e., /u/0o ~ e. Second, we 
limit ourselves to a region near the unitarity limit where 
the boson binding energy is small compared to 0o, i.e., 
£b/0o ~ e so /iB/00 ~ £• The last two terms in £1 give 
tadpoles to the ip and p* fields, which are proportional 
to — 7(70o/co- The condition of cancellation of tadpole 
diagrams will determine the value of the condensate 0o ■ 

Finally, £2 generates two additional vertices for the 
boson propagator as depicted in the last column of Fig. 2, 
which are — illo and —2ifi where 



no(po,p) =Po 



2 ' 



(31) 



These two vertices can be thought of as counterterms to 
avoid double counting of certain types of diagrams which 
are already taken into £0 and £1. For the rest of this 
section, we will consider the unpolarized case with H = 
at zero temperature. 
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B. Power counting rule of e 

We now develop the power counting rule which allows 
us to identify Fcynnian diagrams that contribute to a 
given order in e. Our discussion will extend the power 
counting scheme developed for the unitarity limit in 
Ref. [16] to the region near the unitarity where g'^ /\co\ ^ 

t(f>o. 

Let us first consider Feynman diagrams constructed 
only from Cq and £i, without the vertices from £2- We 
make a prior assumption y(i/0o ~ e, which will be checked 
later, and consider (f)a to be 0{1). Each pair of fcrmion- 
boson vertices is proportional to 5^ ~ e and hence brings 
a factor of e. Also, each insertion of /z ~ e or /zb = 
2/i — .g^/co ^ e brings another factor of e. Therefore, 
the power of e for a given diagram is naively Ng/2 + iV^, 
where Ng is the number of couplings g from C\, and 
iV^ = A^^p + iV^g is the sum of the number of /i insertions 
to the fcrmion line and insertions to the boson line. 

However, this naive power counting does not take into 
account the possibility of inverse powers of e from loop 
integrals. Such a 1/e behavior arises for integrals which 
are ultraviolet divergent at d = 4. Let us identify dia- 
grams which have this divergence. 

Consider a diagram with L loop integrals, Pp fermion 
propagators, and Pb boson propagators. Each momen- 
tum loop integral in the ultraviolet region behaves as 

j dpodp^ j dpspr^p^ (32) 

while each fermion or boson propagator falls at least as 
G{p) ~ p^^ or D{p) ~ or faster. Therefore, the 
maximal degree of divergence that the diagram may have 
is 

r' = 6L-2PF-2PB. (33) 

We call V the superficial degree of divergence. V = 
corresponds to the possibility of logarithmic divergence, 
■D = 2 to quadratic divergence, etc. 

Moreover, following an analysis similar to that done 
in relativistic field theories [25], one finds the following 
relations which involve the number of external fermion 
(boson) lines i^F(B): 

i = (Pf - iV^p) + (Pb - N^^) -Ng + 1, 

(34) 

TV, = (Pf - TV^p) + ^ = 2(Pb - 7V^,) + ^B. 
I 



With the use of these relations, the superficial degree of 
divergence T> can be written as 



V = 6- 2{Ef + Eb + N^), (35) 



which shows that the inverse powers of e appear only in 
diagrams where the total number of external lines and 
chemical potential insertions does not exceed three. This 
is similar to the situation in quantum electrodynamics 
where infinities occur only in electron and photon self- 
energies and the electron-photon triple vertex. 

However, this estimation of V is actually an over- 
estimate: for many diagrams the real degree of diver- 
gence is smaller than given in Eq. (35). To see that, 
we split G{p) into the retarded and advanced parts: 
G(p) = G^ ip) + G'^ip). where G^ (G^) has poles only 
in the lower (upper) half of the complex po plane. It 
is easy to see that the ultraviolet behaviors of different 
components of the propagators are different: 

Gf,{p)^G^,{p)^D''{p)^p-^ (36) 
Gi2(p) ~ G2i(p) ~ p-\ GUp) ~ G^2(P) ~ P"'- 

Note that since the boson propagator D{p) has a pole 
only on the lower half plane of Pq, we have only the re- 
tarded Green's function for the boson D^{p) = 0. 

From these analytic properties of the propagators in 
the ultraviolet region and the vertex structures in C\ 
as well as the relation of Eq. (35), one can show that 
there are only four skeleton diagrams which have the 1/e 
singularity near d = 4. They are one-loop diagrams of 
the boson self-energy [Figs. 3(a) and 3(c)], the tad- 
pole [Fig. 3(e)], and the vacuum diagram [Fig. 3(g)]. We 
shall examine these apparent four exceptions of the naive 
power counting rule of e one by one. 

The first diagram. Fig. 3(a), is the one- loop diagram 
of the boson self-energy. The frequency integral can be 
done explicitly, yielding 



-ina(p) 



dk 



ig' f 

Jk 



1 



{Ek 



G22{k- 



AEk-vEk+ip 



P 
2, 

':){Ek+l 



Ek-^ + Ek+p -po 



£fc+f) , {Ek- 



■MEk, 



-'k+i 



(37) 



El. 



E, 



fc+p +po 



id 
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The integral over k is ultraviolet divergent at d = 4 and 
has a pole at e = 0. Thus H^ip) is 0(1) by itself instead of 
0(e) according to the naive power counting. The residue 
at the pole is 



Po 



f)+0(.), 



(38) 



and is canceled out exactly by adding the vertex Ho in 
£2- Therefore the diagram of the type in Fig. 3(a), when 
combined with the vertex from C2 in Fig. 3(b), conforms 
to the naive e power counting, i.e., is 0(e). 

Similarly, the diagram in Fig. 3(c) representing the bo- 
son self-energy with one /i insertion gives 



dk 



i2n) 



d+l 



Gil (a 



- G 



(39) 




= 0{e) 



(a) 



(b) 




-I- 



+ 




(d) 



(c) 





(e) 



3/2 \ 



(f) 




0(1) 



(g) 



FIG. 3: Four apparent exceptions of the naive power counting 
rule of e [(a), (c), (e), (g)]. The boson self-energy diagram (a) 
or (c) is combined witli the vertex from £,2, (b) or (d), to 
restore the naive e power counting. The tadpole diagram in 
(e) cancels with other tadpole diagrams at the minimum of 
the effective potential. The vacuum diagram (g) is the only 
exception of the naive power counting rule of e, which is 0(1) 
instead of 0(e). 



which also contains a 1/e singularity, and hence is 0(e) 
instead of the naive O(e^). The leading part of this dia- 
gram is 



Jk 4£ 



-2/i + 0(e2), (40) 



which is canceled out exactly by adding the second vertex 
from £2- Then the sum of Figs. 3(c) and 3(d) is again 
O(e^), consistent with the naive power counting. 

The if tadpole diagram with one /x insertion in Fig. 3(e) 

is 



dk 

-M5 / [Gii(fc) - G22{k)] G2i(fc) 



2E^ e0o V 2^ 



(41) 



+ 0(e3/2) 



which is 0(e^/^) instead of the naive 0(e3/2). rpj^ig 
pole diagram should be canceled by other tadpole dia- 
grams of order 0(e^/^) in Fig. 3(f), 



]9_ I mcpo 
' 2 V 27r 



2 (42) 

_,5^ + 0(,3/2)_ 

Co 



The condition of cancellation, Sc + 5f = 0, gives the gap 
equation that determines 0o(m) to the leading order in e. 
The solution to the gap equation is 

2n 2 /2ttV 
00 = — -- — +Oe. 43 
e co\mJ 

When Co < 0, this solution can be written in terms of 
the binding energy as 0o = (2m + eb)/e- The previ- 
ously made assumption fJ./(t>o = 0(e) is indeed correct. 
As we will see later, the cancellation of the tadpole dia- 
grams is automatically achieved by the minimization of 
the effective potential. 

Finally, the one-loop vacuum diagram with one ^ in- 
sertion in Fig. 3(g) also contains the 1/e singularity as 



dk 



{2n) 



-[Gii(fc)-G22(fc)] 



£k 
Ik Ek 



e V 27r 



(44) 



0(e). 



The leading part of this diagram is 0(1) instead of naive 
0(e) and cannot be canceled by any other diagrams. 
Therefore, Fig. 3(g) ~ 0(1) is the only exception of our 
naive power counting rule of e. 

Thus, we can now develop a diagrammatic technique 
for the Fermi gas near the unitarity limit where 5^/|co| ~ 
fi in terms of the systematic expansion over e = 4 — d. 
The power counting rule of e is summarized as follows: 

1. We consider fi/ipo ^ e and regard (/>o as 0(1). 
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0(1) 



FIG. 4: Vacuum diagrams contributing to the effective poten- 
tial up to the next-to-leading order in e. The second diagram 
is 0(1) instead of naive 0(e) because of the 1/e singularity. 



2. For any Green's function, we write down all Feyn- 
man diagrams according to the Feynman rules in 
Fig. 2 using the propagators from £o and the ver- 
tices from Ci. 

3. If there is any subdiagram of the type in Fig. 3(a) 
or 3(c), we add the same Feynman diagram where 
the subdiagram is replaced by the vertex from £2, 
Fig. 3(b) or 3(d). 



4. The power of e for the given Feynman diagram will 
be ©(e^s/^+^f"), where Ng is the number of cou- 
plings g and N^^ is the number of chemical potential 
insertions. 



5. The only exception is the diagram in Fig. 3(g) 
which is 0(1). 

We note that the sum of all tadpole diagrams in Figs. 3(e) 
and 3(f) vanishes if (/)o is the solution of the gap equation. 



C. Effective potential to leading and 
next-to-leading orders 

We now perform calculations to the leading and next- 
to-leading orders in e, employing the Feynman rules and 
the e power counting that have just been developed. To 
find the dependence of 0o on /i, we use the effective po- 
tential method [25] in which this dependence follows from 
the minimization of the effective potential Vcs{4"3)- The 
effective potential at the tree level is given by 



Co 



(45) 



Up to the next-to-leading order, the effective potential re- 
ceives contributions from three vacuum diagrams drawn 
in Fig. 4: fermion loops with and without a /i insertion 
and a fermion loop with one boson exchange. 

The contribution from two one-loop diagrams is 0(1) 
and given by 

^i('^o) = * / jAdTi [inG-i(p) + ^ia^G{p)\ 

(46) 



En 



Using the formula in Eq. (10), we can perform the inte- 
gration over p to find 



Vr 



3) 



(47) 



Substituting d = 4 — e and expanding in terms of e up to 
0(e), this contribution becomes 



= y 



^ ^ 7-3(7 + ln2) ^- 



m(j)o 
2tt 



d/2 



(48) 



^ I l-2(7-ln2) ^- 



where 7 « 0.57722 is the Euler-Maschcroni constant. 

The contribution of the two-loop diagram is 0(e) and 
given by 



^2(^0) =5' 
= 9 



dpdq 



■TY[G{p)<J+G{q)<J^]D{p-q) 



(27r)2d+2 

' ' J^^GlliP)G22iq)D{p~q). (49) 



Performing the integrations over pq and qo , we obtain 

{Ep - ep){Eq - Eg) 



4 Jpq EpEq{Ep + Eq+ £p_q/2) 



(50) 



Since this is a next-to-leading-order diagram and the in- 
tegrals converge at d = 4, we can evaluate the integra- 
tions at c? = 4. Changing the integration variables to 
X — £p/(j)Q, y — £q/(f>o, and cos 6* — p ■ q, the integral can 
be transformed into 



2tt 



X xy sm 



— dx dy de (51) 
Jo Jo Jo 

[f{x) - x][f{y) - y] 

I{x)f{v) [gix, y) - V^cosff] ' 



with/(x) = -I- 1 and g{x, y) = f{x) + f{y) + \{x+y). 
The integration over 9 can be performed analytically to 
lead to 



m4>o 
27r 



d/2 



(52) 



where the constant C is given by a two-dimensional in- 
tegral 



/•OO f^OQ 

C — I dx I dy 
Jo Jo 



[.f{x)~x][J{y)-y] 

f{x)f{y) 



X 9{x, y) - \/g'^{x,y) - xy 
A numerical integration over a: and y gives 
C « 0.14424. 



(53) 



(54) 



Now, gathering up Eqs. (45), (48), and (51), we obtain 
the effective potential up to the next-to-leading order in 
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VMM = Vo{<Po) + Vi(0o) + V2{<Po) 

7-3(7 + ln2) 



Co 



+ 



3 



1 



3Ce } - 



1-2(7 -In 2) 



27r 



d/2 



(55) 



The condition that (j)o minimizes the effective potential 
gives the gap equation; dVes/d4>o = 0. The solution to 
this equation is 



2/Lt 



h = -^[1 + {3C -l + lii2)e\ 



(56) 



2<^, 



e/2 



Co 



d/2 



3C- 



Using the relation of cq with the binding energy of boson 
Eb in Eq. (20), we can rewrite the solution of the gap 
equation in terms of £b up to the next-to-leading order 
in e as 



2/i 



i>o = [1 + {SC - 1 + ln2) e] 



(57) 



£b 
e 



3C 



■ln2 



iln^ii 
2 00 



We note that the leading term in Eq. (57) could be re- 
produced using the mean field approximation, but the 
0(e) corrections are not. The 0{e) correction propor- 
tional to the constant C is a result of the summation of 
fluctuations around the classical solution and is beyond 
the mean field approximation. 



D. Thermodynamic quantities necir unitcirity 

The value of the effective potential VeS at its mini- 
mum determines the pressure P = —Vos{(f)o) at a given 
chemical potential /x and a given binding energy of bo- 
son Substituting the dependence of 0o on /x and £b 
in Eq. (57), we obtain the pressure as 



P = 



1 + 



12 2 



e — 



3£b 



■m(j)o 
2tt 



d/2 



(58) 



The fermion number density N is determined by differ- 
entiating the pressure in terms of ^ as 



N = 



dP 



1 - 27-F21n2 ■ 
1 + e 



2tt 



d/2 



■ (59) 



The nontrivial power of e appears because we have raised 
~ to the power oi 2/d. 

From Eqs. (57) and (60), we can determine the ratio of 
the chemical potential and the Fermi energy fi/sp near 
the unitary limit as 



JJ_ 

£F 



,3/2 



- exp 

_fb_ 
2eF 



;ln< 



8-2e 



1 



3C--(1 



In 2) 



e elne e e^' 
'+2 + ^-2^" 



(61) 



The Ine term in the second line originates from the 
= e^^'^ep into the Ineb/^o term in 



substitution of (pQ 
Eq. (57). We note again that if one is interested in the 
extrapolation to d = 3, one can use Eq. (22) to relate 
£b/£F to akp- The first term in Eq. (61) gives the uni- 
versal parameter of the unitary Fermi gas ^ = m/^fIej^^q' 



= ie3/2 ^ _L,5/2 1^ g _ 0246 6^/2 + . . . . 



(62) 



Here we have used the numerical value for C « 0.14424 
in Eq. (54). The smallness of the coefficient in front of 
e^/^ is a result of a cancellation between the two-loop 
correction and the subleading terms from the expansion 
of the one-loop diagrams around d = i. 

Using Eqs. (58), (59), and (60), the pressure near the 
unitarity limit normalized by epN is given by 



£fA^ d + 2 

Then the energy density E 
from Eqs. (61) and (63) as 



£b 

8£f' 



(63) 



fiN — P can be calculated 



d 



E 



sfN d + 2 



2£f 



e e In e 
4 4 



-In — 

2 £f 



• (64) 



The pressure and energy density in the unitarity limit 
are obtained from ^ via the universal relations depending 
only on the dimensionality of space. 



Then the Fermi energy of a d-dimensional free gas with 
the same number density is 



27r 

£F = — 
m 



1 fd 
2^ (-2 



1 



2/d 



90 



1 - 



l-ln2 



(60) 



E. Quasiparticle spectrum 

The e expansion we have developed is also useful for 
the calculations of physical observables other than the 

thermodynamic quantities. Here we shall look at the 
dispersion relation of fermion quasiparticles. To the 
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■iE = 



FIG. 5: Corrections to the fermion self-energy of order 0(e); 
a fj, insertion and one-loop diagrams. 



leading order in e, the dispersion relation is given by 



LUp (p) = Ep = Y £p + (/"o , which has a minimum at 
zero momentum p = with the energy gap equal to 
A = (jjQ — {2^ + eh)/€- The next-to-leading-order correc- 
tions to the dispersion relation come from three sources: 
from the correction of (/)o in Eq. (57), from a fj, inser- 
tion to the fermion propagator, and from the one-loop 
self-energy diagrams — depicted in Fig. 5. 

Using the Feynman rules, the one-loop diagram of the 
fermion self-energy in Fig. 5 is given by 



/dk 
[a+G{k)a^D{p-k) 



(65) 



+a_G{k)a+D{k ^ p)] 



There are corrections only to the diagonal elements of 
the self-energy and each element is evaluated as 



^iiip) = w' I j^^G,,{k)D{p^k) 

Ek — £fc 

2 Jk Ek{Ek+ek-p/2-po-id)' 



(66) 



and 



dk 



i2n) 



— Gi,{k)D{k-p) 
Ek — £k 

2 JkEk{Ek+ek-p/2+po-i5)' 



(67) 



The dispersion relation of the fermion quasiparticle 
ujf{p) is obtained as a pole of the fermion propagator 
det[G^^{uj,p) + fJ.(J3 — S(ti^,p)] = 0, which reduces to the 
following equation: 



ep + fi~T,ii{uj,p) 

(68) 

To find the 0(e) correction to the dispersion relation, 
we only have to evaluate the self-energy Y,{lo,p) with lj 
given by the leading order solution u) — Ep. Denoting 
Ell (i?p,p) and 'E22{Ep,p) simply by En and E22 and 
solving Eq. (68) in terms of oj, we obtain the dispersion 
relation of the fermion quasiparticle as 




£0 ~ 0(e) 



PL 

2m 



FIG. 6: Illustration of the dispersion relation of the fermion 
quasiparticle in the unitarity limit from the e expansion. 
'^f(p) = \/ (ep — eo)^ + in Eq. (72) is plotted as a func- 
tion of Ep — /2m. Values on both axes are extrapolated 
results to e = 1 in units of the chemical potential /i. 



In order to find the energy gap of the fermion quasi- 
particle, since the minimum of the dispersion relation 
will appear at small momentum Sp ~ /i, we can expand 
T,{Ep,p) around zero momentum p = as J^{Ep,p) — 
E°(0o,O) + E'(0o, 0) £p/0o- Performing the integration 
over fc in E analytically, we have 



and 



Eii(0o,p) =e(2-81n3 + 81n2)(/)o 

+ € I -^ + 81n3-81n2^ £p, 

/,o,p) = e(-2-81n3 + 161n2)(/.o 

e| ---81n3 + 161n2|£ 



(70) 



(71) 



Introducing these expressions into Eq. (69), we find the 
fermion dispersion relation around its minimum has the 
following form: 



ujf{p) ~ a -f- 



(£p — £0 
200 



^ ^ ^(£p-£o)2 + A2. (72) 



Here A is the energy gap of the fermion quasiparticle, 
which is given by 



A = 



^11 + ^2: 



22 



2 (73) 
[1 - (81n3- 121n2)e + 0(e2)] 0o- 



The minimum of the dispersion curve is located at a 
nonzero value of momentum, \p\ = (2to£o)^/^, where 
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Then, introducing the solution of the gap equation (57), 
the energy gap A as a function of the chemical potential 
and the binding energy is given by 

2/x 



A 



[l + (3C-l-81n3 + 131n2)e] 



e 



1+ ( 3C - ^ -81n3 + 131n2- iln^ I e 



= ^(1- 0.345 e) 



(75) 



e 



1 + 0.155 e- 



elne e , A 2u 

h - In 1 + — 

2 2 V £b 



while So is given by 



£0 = 2/x + 



eh 



(76) 



Note the difference with the mean field approximation, 
in which eq = M- 

When So is positive, the fermion dispersion curve has 
its minimum at nonzero value of momentum as in the 
BCS limit, while the minimum is located at zero mo- 
mentum when £o is negative as in the BEC limit. We 
find the former (eq = 2^ > 0) is the case in the unitarity 
limit. The dispersion c;urvc of the fermion quasiparticle 
lof{p) in the unitarity limit £b = is illustrated in Fig 6. 
In particular, the difference between the fermion quasi- 
particle energy at zero momentum and at its minimum 
is given by £q/20o = e/z ~ O(e^). 



F. Location of the splitting point 

Let us consider the situation where the binding energy 
Eb is increasing from zero while the number density N is 
kept fixed. Then A > is held fixed but /x is changing. 
Since to the leading order in e, the chemical potential as 
a function of the energy gap is given by 2/x = cA — e^, 
the location of the minimum of the dispersion curve can 
be written as 

£b 
2 ■ 



£o = eA • 



(77) 



Therefore, Eq decreases as increases. When the bind- 
ing energy reaches eb = 2eA, the minimum of the disper- 
sion curve sits exactly at zero momentum p = 0. This 
point is referred to as a splifMng point [26] . We find the 
splitting point is located on the BEC side of the unitarity 
limit where the binding energy is positive Eb = 2eA > 
and the chemical potential is negative 2^, = — eA. This 
splitting point will play an important role to determine 
the phase structure of the polarized Fermi gas in the uni- 
tary regime as we will study in Sec. IV. 



Although the formalism is based on the smallness of e, we 
see that the next-to-lcading-order corrections are reason- 
ably small compared to the leading terms even at e = 1. 
If we naively use only the leading and next-to-leading- 
order results for ^ in Eq. (62), A in Eq. (75), and sq in 
Eq. (76) in the unitarity limit, their extrapolations to 
e = 1 give for three spatial dimensions 



0.475, 



— w 1.31, 



£^«2. 



(78) 



They are reasonably close to the results of recent Monte 
Carlo simulations, which yield ^ « 0.42, A//i ss 1.2, and 
Eo/m ~ 1-9 [15]. They are also consistent with recent 
experimental measurements of ^, where ^ = 0.51±0.04 [9] 
and ^ = 0.46 ± 0.05 [21]. These agreements can be taken 
as a strong indication that the e expansion is useful even 
at e = 1. 

Our e expansion predicts the behavior of the thermody- 
namic quantities near the unitarity limit. From Eqs. (63), 
(64), and (61), the extrapolations to the three spatial di- 
mensions e = 1 lead to 



2 l£b 
d + 8 £F ' 



E _ d lEb /5 



£F 



In 



£b 
£F 



i^f3-ln£^ 

4 £f V £f 



(79) 



(80) 



(81) 



The pressure, energy density, and chemical potential at 
fixed density are decreasing functions of the binding en- 
ergy Eb with the slope shown above. In Sec. VI, we will 
make a discussion to improve the extrapolation of the se- 
ries expansion by imposing the exact result at = 2 as a 
boundary condition. 

We can also try to determine the location of the split- 
ting point. At this point, the deviation from the unitarity 
point measured by the binding energy per Fermi energy 
is given by 

I = ^[1 + 0(^)1^2 atd = 3. (82) 

Comparing with Eq. (22), one finds afcp w 1 at the split- 
ting point. Since this result is known only to the lead- 
ing order in e, one must be cautious with the numerical 
value. However, certain qualitative features are probably 
correct: the splitting point is located on the BEC side 
of the unitarity limit (a > 0) and the chemical potential 
/Lt/£F « —0.5 is negative at this point. 



IV. PHASE STRUCTURE OF POLARIZED 
FERMI GAS 



G. Extrapolation to e = 1 



A. Effective potential at finite superfluid velocity 



Finally, we discuss the extrapolation of the series ex- 
pansion to the physical case in three spatial dimensions. 



Here we apply the e expansion developed in the previ- 
ous section to the unitary Fermi gas with unequal den- 
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sities for the two fcrmion components. Wc put special 
emphasis on investigating the phase structure of polar- 
ized Fermi gas in the unitary regime, which has a di- 
rect relation with the recent measurements in atomic 
traps [20-23]. For this purpose, we generalize our for- 
malism to take into account the possibility of a spa- 
tially varying condensate where {4>{x)) = e^^'^^'^^o with 
Vs being the superfluid velocity. The factor e^'™''^''" in 
front of (fiQ in the Lagrangian density can be absorbed 
by making the Galilean transformation on the fermion 
field as ipcr{x) — » e''^'"''''^tjj„{x) and the boson field as 



"(fiix). Accordingly, the fermion propaga- 



tor in Eq. (29) is modified as 
G{po,p) = 



{PO + H -p- Vs)'^ - {Sp + SmvJ'^ - <Po + 

Po+ H + 

Sp—mVs 00 
— 00 Po + H — Sp+mv, 

(83) 

while the boson propagator in Eq. (30) is modified as 



D{po,p)={po-'-^^ + i6y 



(84) 



The value of the superfluid velocity Vg is determined 

by minimizing the Vs dependent part of the effective po- 
tential Vh{vs)- As far as the polarization H is sufficiently 
small compared to the energy gap A of the fermion quasi- 
particle, Vh{Vs) to the leading order in e is given by the 
fermion one-loop diagram with the propagator in Eq. (83) 
as follows: 



+ 



Po- 



(85) 



Since it will turn out that Smv^ is 0(e^^), we can expand 
the integrand in terms of emv^/<j>o to lead to 



27r 



(86) 



Using the result on the fermion number density A'' in 
Eq. (59), this part can be rewritten as Vh</!\{vs) = 
Nsmve ' which represents the energy cost due to the pres- 
ence of the superfluid flow. 

li H —p Vs reaches the bottom of the fcrmion quasipar- 
ticle spectrum A + {£p + £mv^—So)'^ /"^(l^o, ^//(^s) receives 
an additional contribution from the filled fermion quasi- 
particles. Since Smvs ^ is small, we can neglect it in 
the quasiparticle energy. Then the Vs dependent part of 
the effective potential is given by 



P Vs 



wf(p) 



(87) 



where wc have introduced a notation {x)y = 9{x) and 
ujf{p) is the fermion quasiparticle spectrum ojf{p) = A-|- 
(erp— £o)^/2</>o derived in Eq. (72). The effective potential 
Vh{vs) to the leading order in e has the same form as that 
studied based on the effective field theory in [26]. 



B. Critical polarizations 

Now we evaluate the effective potential Vh{vs) at d = 4 
as a function of Vg and H. Changing the integration 
variables to z = £p/<po and w = p ■ v^, we have 



m(f>o 

277 



— / dz dw z\/ 1 — w'^ 

Jo J^l 



(z — Zn)'^ I 

X I - A - — ;^-^</>o - z4>Qejnv,W 



(88) 



where = £o/?!'o ~ £• We can approximate ^Jz^qe^^ 
by \/ ZQ^t^Smv^ because the difference will be O(e^) and 
negligible compared to itself \J ZQ<^Q£mv^ ~ ^ ■ Then the 
integration over z leads to 



27r 



e<po 



^ d'w\f 



32 
iJ- A 



-zo 



2-\/ ^mvAo^Ki 



3/4 



3/2' 



(89) 



By introducing the dimensionless variables x and h as 



= X I — e I z, 



and 



7 

' 



'0 5 



(90) 



(91) 



we can rewrite the effective potential Vff(t's) = V/i(a;) in 
the simple form as 







(92) 



- a;3/2 y ' dw v/r^ Q - 



w 



Z/2 



Now we see that if the superfluid velocity exists, Emv^ ~ 
e^^ and the Vs dependent part of the effective potential 
isO(eiO). 

Numerical studies on the effective potential Vh{x) as a 
function of x show that there exists a region of h, h\ < 
h < h2, where Vh{x) has its minimum at finite x. These 
two critical values are numerically given by 



hi = -0.00366 and /12 = 0.0275. 



(93) 



Correspondingly, we obtain the critical polarizations nor- 
malized by the energy gap as 



Hi 



1 = 1-0.0843.^(1) , 



(94) 



13 



and 



^ = 1 + 0.634 
A VA 



(95) 



Here we have replaced 0o by A because they only differ 
by 0(e) [Eq. (73)]. The region for the supcrfluid phase 
with the spatially varying condensate is H2 — Hi ~ e^, 
where the superfluid velocity Vg is finite at the ground 
state. 

As the polarization increases further H > H2, the 
superfluid velocity disappears. If iJ < wf(0), fermion 
quasiparticles which have momentum ljf{p) < H are 
filled and hence there exist two Fermi surfaces, while 
there is only one Fermi surface for H > uJpiO) (see Fig. 6). 
Therefore, from the quasiparticle spectrum derived in 
Eq. (72), the polarization for the disappearance of the 
inner Fermi surface H3 is given by 



H3 
A 



wf(0) 



A 2Va, 



(96) 



Here the location of minimum in the fermion quasipar- 
ticle spectrum eq is related with the binding energy Sh 
near the unitarity limit via Sq = eA — eb/2 [Eq. (77)]. 
The critical polarizations Hi/ A and H2/A as functions 
of Sh/eA are illustrated in Fig. 7. 



C. Phase transition to normal Fermi gas 



The polarization dependent part of the pressure Ph ~ 
is small and negligible in the region considered here 
- A ~ e. 

The pressure of the superfluid state should be com- 
pared to that of the normal state with the same chemical 
potentials. Since the phase transition to the normal state 
happens at if ~ 00 ^ ~ e^o, only one component of 
fermions exists in the normal state. Therefore, the in- 
teraction is completely suppressed in the normal Fermi 
gas and its pressure is simply given by that of a free 
Fermi gas with a single component as follows: 



-Pn = / (MT -Sp)> = 



{H + n)i+^ 
r(f + 2) V27r 



(100) 



The phase transition of the superfluid state to the nor- 
mal state occurs at H = H^ where the two pressures 
coincide P = P^. From Eqs. (99) and (100), the critical 
polarization He satisfies the following equation: 



Hr = 



Eb „ 2 + ln2 ■ 
400 6 



1 + ^-Ce- ?±i^i^e + (81n3 - 121n2)e 
4A 6 ^ ^ 



(101) 
A, 



where we have substituted the relation between the con- 
densate 4>Q and the energy gap A at zero polarization in 
Eq. (73). Defining a number a 0(1) by 



Next we turn to the phase transition to the normal 
Fermi gas which occurs at i7 — A ^ e. Since the region 
for the superfluid phase with the spatially varying con- 
densate is ifi,2 — A ~ e^, we can neglect the superfluid 
velocity here. We can also neglect the possibility to 
have two Fermi surfaces where H3 — A ^ e^. Then the 
contribution of the polarized quasiparticles to the effec- 
tive potential is evaluated as 



Vh{0) 



[H-u;f{p)V 



{H 



A)l 



200 



"^0C 
277 



d/2 



(97) 

where we have neglected the higher-order corrections due 
to the shift of the location of minimum in the dispersion 
relation £0 ~ e- According to the modification of the ef- 
fective potential in Eq. (55) to Vcfi{(t>Q) + Vh , the solution 
of the gap equation in Eq. (57) becomes 0o ^ 0o + 4>h, 
where 



6^ = _(iJ_A)> + 



{H-A)l 
200 



(98) 



Then the pressure of the polarized Fermi gas in the su- 
perfluid state is given by P = — 14ff(0o + 0ff) — V^r, which 
results in 



1+('17_3^_7+In2\ 3£b 



2-K 



^0{e{H-A)l, {H-A)l- 



dji 



(99) 



a = C^ 2+^ln2 _ (glnS - 121n2) = 0.122, (102) 

the critical polarization normalized by the energy gap at 
zero polarization is written as 



^= 1 I 2 



(103) 



If the binding energy is large enough Sb/eA > 4cr = 
0.488, the superfluid state remains above H = A. The 
superfluid state at A < iJ < iJc is referred to as a gap- 
less superfluid state. The fermion number densities for 
two different components are asymmetric there and the 
fermionic excitation does not have the energy gap. 

In particular, at the unitarity limit where Sh/A = 0, 
the critical polarization is given by 



He 
A 



= l-ecr = 1- 0.122 e. 



UL 



At the splitting point where Sb/A = 2e, it is 



-He 
A 



1 



SP 



0.378 e. 



(104) 



(105) 



The extrapolations to three spatial dimensions e = 1 
give the critical polarizations at the two typical points 
as Fc/A|uL = 0.878 and Fc/A|sp = 1-378. At unitar- 
ity, hence, there is no gaplcss superfluid phase due to the 
competition with the normal phase. On the other hand. 
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FIG. 7: Schematic phase diagram of the polarized Fermi gas 
near the unitarity limit from the e expansion. The phase dia- 
gram can be divided into four phases, I: the gapped superfluid 
phase, II: the polarized normal phase. III: the gapless super- 
fluid phase, and IV: the phase with the spatially varying con- 
densate. The inset is the magnification of the region around 
the splitting point (SP). The phase IV appears in the narrow 
region represented by the thick line between the phases I and 
III. 

near the splitting point, the normal phase is not com- 
petitive compared to the gapless phases and the gapless 
superfluid phases stably exist there. The phase bound- 
ary between the superfluid and normal phases He/ A as 
a function of Sb/eA is illustrated in Fig. 7. 

D. Phase structure near the unitarity limit 

The schematic phase diagram of the polarized Fermi 

gas in the unitary regime is shown in Fig. 7 in the plane 
of H and eb/e for the fixed energy gap A at zero polar- 
ization. The critical polarization He/ A divides the phase 
diagram into two regions; the superfluid phase at H < He 
(I and III) and the normal phase sX H > He (II). The 
Fermi gas in the normal phase is fully polarized near the 
unitarity limit because H ^ ^. The phase transition of 
the superfluid state to the normal state is of the first or- 
der, because of the discontinuity in the fcrmion number 
density which is 0{e~^) in the superfluid phase while it 
is 0(1) in the normal phase. 

The superfluid phase can be divided further into two 
regions; the gapped superfluid phase at iJ < A (I) and 
the gapless superfluid phase at A < H < He (HI). On the 
BEG side of the splitting point (SP) where £b/eA > 2, 
the phase transition from the gapped phase to the gap- 
less phase is of the second order because a disconti- 
nuity appears in the second derivative of the pressure 
d'^P/dH'^ ~ e0{H - A) [see Eq. (99)]. On the other 
hand, on the BCS side of the splitting point where 
4cr < £b/eA < 2, there exists the superfluid phase with 
the spatially varying condensate (IV) at Hi < H < H^ 
between the gapped and gapless phases. This phase ap- 
pears only in the narrow region where H2 — Hi ~ . The 



eH/ei, 

* H = H, 




1 1 1 1 >■ — 

-0.5 0.5 1 £b 

FIG. 8: Schematic phase diagram of the polarized Fermi gas 
near the unitarity limit in the H-^ plane for the fixed binding 
energy > 0. The phases I, II, and III are the same as in 
Fig. 7. The horizontal dashed line from right to left tracks 
the chemical potentials {n, H) in trapped Fermi gases from 
the center to the peripheral. 



phase transitions at H = Hi, H2 are of the flrst order. 

In actual experiments using fermionic atoms, unitary 
Fermi gases are trapped in optical potentials V{r) [20- 
23]. In such cases, the polarization H and the scattering 
length or Cb are constant over the space, while the effec- 
tive chemical potential iJ,—V{r) decreases from the center 
to the peripheral of the trapping potential. For the pur- 
pose of comparison of our results with the experiments 
on the polarized Fermi gases, it is convenient to visualize 
the phase diagram in the plane of H and fi for the fixed 

In Fig. 8, eHe/sh and eA/s^ are plotted as functions 
of yu/eb for the fixed binding energy Sh > 0- The su- 
perfluid phase and the normal phase are separated by 
the line of the critical polarization H = H^ while the 
gapped and gapless superfluid phases are separated by 
the line of the energy gap at zero polarization H = A. 
The intersection of the two lines. He = A, is located 
at (/x/£b,eiJ/£b) = ((1 -4a)/8a, 1/4(7) = (0.524,2.05), 
while the splitting point in the phase diagram is at 
(/i/£b, eff/£b) = (^1/4,1/2). The horizontal dashed line 
from right to left tracks the chemical potential /x in the 
trapped Fermi gas from the center to the peripheral. If 
the polarization is small enough compared to the binding 
energy eH/sh < 2.05, there exists the gapless superfluid 
phase (III) between the gapped superfluid phase (I) and 
the normal phase (II). When the polarization is above 
the splitting point 0.5 < eH/sh < 2.05, the phase with 
the spatially varying condensate (thick lino between I and 
III) will appear between the gapped and gapless super- 
fluid phases. 

Assuming that the above picture remains qualitatively 
valid in three spatial dimensions, we thus found that the 
physics around the splitting point is the same as argued in 
Ref. [26]. However, due to the competition with the nor- 
mal phase, the gapless phases disappear at some point, 
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probably before the unitarity is reached. The point where 
this happens can be estimated from our calculations to 
be 



iG 



= — ^ 



2a = 0.494. 



(106) 



The gapless phase with the spatially varying condensate 
exists in a finite range 0.494 < (afcp)^^ < 1. The hy- 
pothesis of Ref. [26] that the region of the gapless phase 
with the spatially varying condensate is connected with 
the FFLO region on the BCS side is not realized in the 
e expansion. 



V. EXPANSION AROUND TWO SPATIAL 
DIMENSIONS 

A. Lagrangian and Feynman rules 

In this section, we formulate the systematic expansion 
for the unitary Fermi gas around two spatial dimensions 
in a similar way as we have done for the e — A — d expan- 
sion. Here we start with the Lagrangian density given 
in Eq. (24) limited to the unpolarized Fermi gas in the 
unitarity limit where H — and I/cq = as follows: 



(107) 



Then we expand the field (j) around its vacuum expecta- 
tion value (/<o as 



I'D + 9V, 



27re 
m 



1/2 - 

(^) ' 



where the effective coupling g ~ e^/^ in Eq. (18) was 
introduced. The extra factor (m/^/27r) '^^'^ was chosen so 
that the product of fields ip* has the same dimension as 
the Lagrangian density. 

Then we rewrite the Lagrangian density as a sum of 
three parts, C = Co + Ci + C2, where 

Co = (^idt + + H(T3 + a+(j}o + * , 

(109) 

Ci^ -ip*ip'^g^'^a+^Lp + g^^a-^Lp*, (110) 

C2 = ip*ip. (Ill) 

The part Cq describes the gapped fermion quasiparticle, 
whose propagator is given by 



G{po,p) 



1 



p2 _ + i5 



t>o Po 



90 

Sp + p, 



(112) 



iga. 



FIG. 9: Feynman rules for the expansion around two spatial 
dimensions from the Lagrangian density in Eq. (109). The 
first line gives propagators, while the second line gives ver- 
tices. 



with Ep = yj {Sp — + 4>Q being the usual gapped 
quasiparticle spectrum in the BCS theory. 

The second part Ci describes the interaction between 
fermions induced by the auxiliary field Lp. The first term 
in Ci gives the propagator of the auxiliary field Lp, 



D{po,p) 



-1, 



(113) 



and the last two terms give vertices coupling two fermions 
with ip. If we did not have the part £2, we could integrate 
out the auxiliary fields ip and ip* to lead to 

Ci «'ta_^' = gVjV'IV'iV'T: (114) 

which gives the contact interaction of fermions with the 
small coupling ~ e as depicted in Fig. 1. The vertex 
in the third part C2 plays a role of a counterterm so as 
to avoid double counting of a certain type of diagrams 
which is already taken into Ci as we will see below. The 
Feynman rules corresponding to these Lagrangian densi- 
ties are summarized in Fig. 9. 



B. Power counting rule of e 

We can construct the similar power counting rule of e 
as in the case of the expansion around four spatial di- 
mensions. Let us first consider Feynman diagrams con- 
structed only from Cq and Ci, without the vertices from 
£2- We make a prior assumption 0o//^ e~^^'^, which 
will be checked later, and consider p, to be 0(1). Since 
e"^/*^ is exponentially small compared to any powers of 
e, we can neglect the contribution of 0o when we expand 
physical observables in powers of e. Then, since each pair 
of fermion and ip vertices brings a factor of e, the power 
of e for a given diagram is naively Ng/2, where Ng is 
the number of couplings g from Ci. However, this naive 
power counting does not take into account the possibility 
of inverse powers of e from loop integrals which are loga- 
rithmically divergent at d = 2. Let us identify diagrams 
which have this divergence. 

Consider a diagram with L loop integrals, Pp fermion 
propagators, and Pb internal auxiliary field lines. Each 
momentum loop integral in the ultraviolet region behaves 
as 



j dpodp^ J 



dpsp ~ p'^ 



(115) 
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while each fermion propagator falls at least as G{p) ~ 

or faster. Therefore, the superficial degree of divergence 

the diagram has is 

V^4:L-2Pf. (116) 

Using the similar relations to Eq. (34), which involve the 
number of external fermion (auxiliary field) lines i?F(B) , 



L^Pf + Pb-N. + I, 



Na=PF 



Ef 



= 2Pb + Eb, 



(117) 



the superficial degree of divergence can be written as 

V = 4:- Ef- 2Eb. (118) 

This shows that the inverse powers of e appear only in 
diagrams where Ep + 2Eb < 4 is satisfied. Moreover, 
from the analytic properties of G{p) in the ultraviolet 
region discussed in Eq. (36), one can show that there 
are only two skeleton diagrams which have the 1/e sin- 
gularity near d — 2. They are one-loop diagrams of the 
boson self-energy [Fig. 10(a)], and the ip tadpole diagram 
[Fig. 10(c)]. We shall examine these apparent two excep- 
tions of the naive power counting rule of e. 

The boson self-energy diagram in Fig. 10(a) is evalu- 
ated as 



-9 



dk 



{2Tr)d+ 



^-1 



G 



22 



k-P 



k 2ek - {po 



2^P 



+ 2fi + iS) ' 



(119) 



where we have neglected the contribution of (f>Q. The 
integral over k is logarithmically divergent at d = 2 and 
has a pole at e = 0. Thus fii,{p) is 0(1) instead of 0(e) 
according to the naive power counting. The residue at 
the pole can be computed as 

i(p) = -3^ / 7^ : r — — FT + ■ ■ ■ 

Jk 2ek - [Pa - + 2^ + id) 

= l + 0(e-), (120) 



^3 




FIG. 11: The tadpole diagram which gives the medium-effect 
correction to the gap equation. 



and is canceled out exactly by adding the vertex Hq — ~1 
in £,2- Therefore the diagram of the type in Fig. 10(a) 
should be combined with the vertex from £2 in Fig. 10(b) 
to restore the naive e power counting result, i.e., 0(e). 

Similarly, the tadpole diagram in Fig. 10(c) also con- 
tains a 1/e singularity. The requirement that this tadpole 
diagram should vanish by itself gives the gap equation to 
the leading order that determines the condensate 0o as 
we will see in the succeeding subsection. 



C. Gap equation 

The condensate (po as a function of the chemical poten- 
tial n is determined by the gap equation which is obtained 
by the condition of the disappearance of all tadpole dia- 
grams. The leading contribution to the gap equation is 
the one-loop diagram drawn in Fig. 10(c), which is given 
by 



dk 



^1 (2^)rf+i G2i{k)^t9 



2Ek 



By changing the integration variable to z 
obtain 



(121) 
rfc//i, we 



dz 



^e/2 



^#0 (^) ' 

2m r(f) Jo ^{z-ir + id^o/f^y 



(122) 



The integration over z in the dimensional regularization 
can be performed to lead to 




O(e-) 



(b) 




(c) 



FIG. 10: Two apparent exceptions of naive power counting 
rule of e, [(a), (c)]. The boson self-energy diagram (a) is com- 
bined with the vertex from £.2 (b) to restore the naive e power 
counting. The condition of disappearance of the tadpole di- 
agram (c) gives the gap equation to determine the value of 
condensate (j)o- 



^'2tF) 



lnH/^-i + 0(e-) 
00 e 



(123) 



The first term ln2/x/(/)o originates from the singularity 
around the Fermi surface as is well known as the Cooper 
instability, while the second term 1/e is from the loga- 
rithmic singularity of the fc integration at d = 2 in Si. 
Solving the gap equation Si = 0, we obtain the conden- 
sate as (/)o = 2iie~^l^ . Note that this result is equivalent 
to that obtained by the mean field BCS theory. 

It is known that the preexponential factor in the mean 
field result 0o = 2/ie^^/' is modified due to the ef- 
fects of medium [27, 28]. In the language of the tad- 
pole diagrams, the corresponding modification to the gap 
equation comes from the three-loop diagram S3 depicted 
in Fig. 11. The diagram, which seems proportional to 

~ e^, turns out to give the 0(1) correction to the gap 
equation. 
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Using the Fcynman rules, the tadpole diagram in 
Fig. 11 is given by 



dk dp dq 



X Giiip-k)G2i{q)G22{q + k). 
Since the second term in the product, 

1 <l>^ 



(124) 



Gn(p) G22(p) 



(125) 



gives only the 0(e) correction to tlie gap equation, we 
can neglect it for the current purpose. Then we perform 
the integrations over po and g'o to result in 



"3 



/dki 



dpdq 



f'^+'^ AEpEq 

Ep — ko+ £p-k — M 



(126) 



+ 



Ep + £r 



{Ep-kor-E^_^ {ko + Ep^-E^^^ 



+ 



E„ 



ko+ Eg- Eq+k + /i 

[ko + Egf-E^^^ ' {ko-Eg)^-E^_^ 



Because of the factor 1/EpEq in the integrand, the inte- 
grations over p and q are dominated around the Fermi 
surface where ep(q) ~ /x and hence -Ep(q) ~ (/)o <C /x. 
Keeping only the dominant part in the integrand, we can 
write the integral as 



^3 - g" 



dk dp dq 00 
(27r)3rf+i 4^ 
1 



(127) 



(fco + Sp-k - m)(^0 + £q+fe - M) 



Now the integration over ko can be performed easily to 

lead to 



■=■3 = 



pq ^p^q 



Sq+k Sp—k 



(128) 



If we evaluated the fe integration in the second line at 

d = 3, we would obtain the static Lindhard function 
representing the medium-induced interaction. 

Here we shall evaluate S3 at d = 2. Changing the 
integration variables to Sp, Sq, Sk, cosXp = k p, cosXg = 
k q, and performing the integrations over Sp and Sq, we 
obtain 



S3--^^0o(£)^^21n^ 



2 VZTT/ \ (pQ ^ 

" dsk dxp dxg 0{£q+k - IJ) 0{n - Ep-k) 



(129) 



Sq+k 



The range of the integrations over Xp a-nd Xq are from 
to TT. Since 0o//x e~^/'^, we have (21n/i/(/)o)^ ~ (2/e)^ 
which cancels coming from the four vertex couplings 
Finally the remaining integrations can be performed 
as follows: 



dskdxpdxq 0{eq+k - IJ')Oin- Ep-k) 



r dskdx 

J TT- 



Sq+k Sp—k 



Jo Jo Jo J4ticos^Xq Jo ^ A/2 7^ 



2 VMEfe (cos Xp + COS Xq) 



(130) 



which gives S3 as 



receives the 0(1) correction due to S3, is modified as 



S3 ~ -ig(f)o 



m 
2^' 



(131) 



Consequently, the gap equation ^1+^3 =0, which 



ln^-i-l + O(e) = 0. 
90 e 



(132) 
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0(1) 




FIG. 13: One-loop diagrams contributing to the fermion self- 
energy to order 0{e). 



FIG. 12: Vacuum diagrams contributing to the pressure up 
to the next-to-leading order in e. 



The solution of the gap equation becomes 



2/z exp 



1 



1 + 0{e] 



[l + 0(e-)]e- 



-l/e 



(133) 

where the value of condensate is reduced by the factor 
e « 2.71828. The reduction of the preexponential factor 
due to the medium effects is known as the Gor'kov-Melik- 
Barkhudarov correction at d = 3 theories [27, 28]. 



D. Thermodynamic quantities 

The value of the effective potential VcS at its mini- 
mum determines the pressure P = — t4ff('^o) ^ given 
chemical potential /i. Since the energy gain due to the 
superfluidity <^g ~ e"^/*^ is exponentially small compared 
to any power series of e, we can simply neglect the con- 
tribution of 4>o to the pressure. Up to the next-to-leading 
order, the effective potential receives contributions from 
two vacuum diagrams drawn in Fig. 12: fermion loops 
without and with an exchange of the auxiliary field. The 
one-loop diagram is 0(1) and given by 



yi(0) = -2 / (Ai-£p)> = - 



2m(^)' _ 



r(| + 2) 



free ; 



(134) 

which represents the contribution of free fermions to the 
pressure. The two-loop diagram is 0(e), which represents 
the density-density correlation as 



^2(0) = -f 



dpdq 
(27r)2'i+2 

J V 



Gn{p)G22{q) 



(135) 



Thus we obtain the pressure up to the next-to-leading 
order in e as 



P={l + e)Pi 



free ■ 



(136) 



Accordingly, the fermion number density is given by 
N — dP/dfi = (1 -I- e) A^froc- The Fermi energy is ob- 
tained from the thermodynamics of free gas in d spatial 
dimension as 



2n 

ep = — 

TO 



2 \ 2 



1 N 



2/d 



(1 



(137) 



which yields the universal parameter of the unitary Fermi 
gas from the e expansion as 



C=-^-l-e + 0(e'). 



(138) 



The 0{e^) correction to f was recently computed to ob- 
tain 1.52 [29]. 



E. Quasiparticle spectrum 

To the leading order in e, the dispersion relation of 
the fermion quasiparticle is given by wf(p) = Ep = 
\/ {sp — ^)2 -|- , which has the same form as that in the 
mean field BCS theory. There exist the next-to-leading- 
order corrections to the fermion quasiparticle spectrum 
from the one-loop self-energy diagrams —iT,{p) depicted 
in Fig. 13. These corrections are only to the diagonal 
elements of the self-energy and each element is evaluated 
as 



dk 



Sii(p) = -i.9 y(^^^G22(fc) 
= -9^ / ^{fJ- - £fe) = -e^, 

Jk 



(139) 



and 



dk 



6{n ~ Ek) = e^- 



(140) 



To this order, the self-energy is momentum independent, 
which effectively shifts the chemical potential due to the 
interaction with the other component of fermions. 

By solving the equation det[G'~^(ijj,p) — S] = in 
terms of w, the dispersion relation of the fermion quasi- 
particle up to the next-to-leading order is given by 



^f{p) = \/(£p - M - ^t-Y 



(141) 



The minimum of the dispersion curve is located at a 
nonzero value of momentum, |p| = (2to£o)^^^, where 



£o = (1 + e)M- 



(142) 



The location of the minimum coincides with the Fermi 
energy in Eq. (137), = ep, in agreement with the Lut- 
tinger theorem [30]. The energy gap A of the fermion 
quasiparticle is given by the condensate, 



00 = — e / 
e 



(143) 
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F. Extrapolation to e=l 



Now we discuss the extrapolation of the series expan- 
sion over e — d — 2 to the physical case in three spatial 
dimensions. In contradiction to the case of the e = A — d 
expansion, the coefficients of the 0(e) corrections are not 
small compared to the leading terms. If we naively ex- 
trapolate the leading and next-to-leading-order results 
for ^ in Eq. (138), A in Eq. (143), and Eq in Eq. (142) to 
e = 1, we would have 





2 



^ 4 ii-i 



0, —« 0.271, 



(144) 



which are not as good as the extrapolations in the ex- 
pansions over e in Eq. (78). This may be related to the 
fact that the series expansion over e = d — 2 is an asymp- 
totic series and is not Borel summable. Thus, instead 
of naively extrapolating the e expansions to d = 3, we 
use them as boundary conditions to improve the series 
summation of the e expansions in Sec. VI. 



G. Large-order behavior 



FIG. 14: An nth order diagram at d = 2 which contributes to 
the pressure as n! by itself. The countervertex — jflo = i for 
each bubble diagram is understood implicitly. 



can be evaluated as 



1+5^ 



= — e In 



(gfc+E - m) - ^(m - gfe-f) 
k 2ek - [pq - \ep -1-2^-1- i5) 




(146) 



Here we show that the series expansion over e = d — 2 
is not convergent, but an asymptotic series. Its physical 
reason is obvious: The unitary Fermi gas at d < 2 is just a 
free gas while the ground state at d > 2 is the superfluid. 
Therefore, the radius of convergence in the expansion 
around two spatial dimensions is zero. In the language 
of the perturbation theory developed in this section, the 
asymptotic nature of the expansion over e — d~2 can be 
understood by the existence of a type of diagrams which 
grows as nl by itself at order e". Such a factorial con- 
tribution originates from the low momentum integration 
region which resembles the infrared renormalon in the 
relativistic field theories [31, 32]. 

The n + 1-loop diagram contributing to the effective 
potential as nl at 0(e") is depicted in Fig. 14, which is 
written as 



Vn 



+ 1 



dp 



(27r)'^+i 

6'(£fc+p - m) - e{fl - Ek-R, 



(145) 



k '2ek - {po 



Is 



+ 2ii + iS) 



where feg integrations in each bubble diagram are already 
performed. Note that +1 in the bracket comes from the 
countervertex ^Hq ~ 1. Since the n\ contribution comes 
from the infrared region of p, we concentrate on the in- 
tegration region where |po| ^ M ^-nd Sp <C /U. Rotating 
the contour of the po integration to the imaginary axis, 
Po —^ ipo (Wick rotation) , the k integration in the bracket 



Here • • • includes higher order terms in e and regular 
dependences on pQ and Sp. The logarithmic dependence 
on Po and Sp appears due to the singularity of the k 
integration near the Fermi surface £fe ~ /i in the limit p — > 
0. Such a logarithmic term dominates the p integration 
when p is small. Collecting the logarithmic term from 
each bracket, the po and p integrals in Eq. (145) at their 
low momentum regions contribute to Vn+i as 



dpodp 



-e In 




(I) 



3) r(n), 



which behaves as the factorial of n. 



(147) 



Because the large-order contribution grows up with the 
same sign, it produces the singularity on the real positive 
axis of the Borel transform. Thus the expansion over e = 
d — 2 is not Borel summable and its Borel transformation 
will have the inevitable ambiguity '-^ e"'^/'^, indicating 
that the ground state where the expansion is performed 
is unstable. Finally, it is important to note that the 
discussion here is not applicable to the expansion around 
d = 4 because the singularity near the Fermi surface 
is cured by the existence of the large condensate 4>o ^ 
fi. In the Appendix, we explicitly demonstrate that the 
diagrams in Fig. 14 do not yield a factorial growth in the 
coefficients of the e expansion. 
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VI. MATCHING OF EXPANSIONS AROUND 
FOUR AND TWO SPATIAL DIMENSIONS 

As we have mentioned previously, we shall match the 
two expansions around d = 4 and d = 2, which were stud- 
ied in Sees. Ill and V, respectively, in order to extract 
results at d = 3. We use the results around two spa- 
tial dimensions as boundary conditions which should be 
satisfied by the series summations of the expansions over 
e = 4 — d. Because we do not yet have a precise knowl- 
edge of the large-order behavior of the expansion around 
four spatial dimensions, we assume its Borel summability 
and employ Fade approximants. 

Let us demonstrate the matching of the two expan- 
sions by taking ^ as an example. In Ref. [17], the linear 
interpolation between the exact values at d = 2 = 1) 
and d = 4 = 0) was discussed to yield — 0.5 at 
d = 3. Now we have series expansions around these two 
exact limits. The expansion of ^ in terms of e = 4 — d is 
obtained in Eq. (62). Assuming the Borel summability 
of the e expansion, we write ^ as a function of e in the 
form of the Borel transformation, 



^(e) = — exp(^ 



Jo 



(148) 



where we factorized out the trivial nonanalytic depen- 
dences on e explicitly. (t) is the Borel transform of the 
power series in ^(e), whose Taylor coefhcients at origin 
are given from the e expansion of ^ as 



B^{t) = l-[3C--il 



ln2) U-l- 



(149) 



In order to perform the integration over t in Eq. (148), 
the analytic continuation of the Borel transform B^ (t) to 
the real positive axis of t is necessary. Here we employ 
the Fade approximant, where B(^{t) is replaced by the 
following rational functions: 



Bdt) 



l+Pit 



- Put 



M 



l + qit^ + 



(150) 



From Eq. (149), we require that the Fade approximants 
satisfy pi — qi — — 3C -I- |(1 — In 2). Furthermore, we 
incorporate the results around two spatial dimensions in 
Eq. (138) by imposing 



1 



(151) 



on the Fade approximants as a boundary condition. 
Since we have three known coefficients from the two 
expansions, the Fade approximants [M/N] satisfying 
M -I- TV = 3 are possible. Since we could not find a solu- 
tion satisfying the boundary condition ^(2 — e) — 1 — e for 
[M/N] = [2/1], we adopt three other Fade approximants 
with [M/N] = [3/0], [1/2], [0/3], whose coefficients p„, 
and qn are determined uniquely by the above conditions. 

Figure 15 shows the universal parameter ^ as a function 
of the spatial dimension d. The middle three curves show 




FIG. 15: (Color online) The universal parameter ^ as a func- 
tion of the spatial dimensions d. The upper solid curve is from 
the expansion around d = 4 in Eq. (62), while the lower solid 
line is from the expansion around d = 2 in Eq. (138). The 
middle three curves show the different Borel-Pade approxima- 
tions connecting the two expansions. The diamond indicates 
the result ^ ~ 0.42 from the Monte Carlo simulation [15]. 



^ in the different Fade approximants connecting the two 
expansions around d = 4 and d = 2. These Borel-Fade 
approximations give ^ = 0.391, 0.364, and 0.378 at d = 3, 
which are small compared to the naive extrapolation of 
the e expansion to d = 3 (^ ^ 0.475). 

The Fade approximant above, however, almost cer- 
tainly needs serious modification, since we can show that 
in the expansion of ^(e), there exist nonanalytic terms at 
sufficiently higher orders in e. 



g3/2 

C(e) = — exp 



elne 



8-2e^ 

(1 - 0.0492e + # £2 + # 1^ g _ 



(152) 



The nonanalytic term at 0{e^) compared to the lead- 
ing order originates from a resummation of boson self- 
energies to avoid infrared singularities. An understand- 
ing of the analytic structure of high-order terms in the 
perturbation theory around d = 4 is currently lacking. 



VII. SUMMARY AND CONCLUDING 
REMARKS 

We have developed the systematic expansion for the 
Fermi gas near the unitarity limit treating the dimen- 
sionality of space d as close to four. To the leading and 
next-to-leading orders in the expansion over e = 4 — d, 
the thermodynamic quantities and the fermion quasipar- 
ticle spectrum were calculated as functions of the binding 
energy of the two-body state. Results for the physical 
case of three spatial dimensions were obtained by extrap- 
olating the series expansions to e = 1. The pressure, en- 
ergy density, and chemical potential in the unitary regime 
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at fixed density were found to be 

P 2 ^ e £b 



epN d + 4 2£F 

E d ifb /, e elne e Sh 
■ ? - 7, — 1 + 7 "I ■ ~ 1^ — 



(153) 



£fN d + 2eF \ 4 4 2 ep 

(154) 

where 

3/2 5/2 

C = — + — lne-0.0246e^/2^O(e''/2) « 0.475 (156) 

is the universal parameter of the unitary Fermi gas. 
The fermion quasiparticle spectrum around its minimum 
in the unitarity limit was given by the form ujf{p) — 
\/ {sp — so)^ + with the energy gap 

A 2 

_ = - -0.691 + 0(e) « 1.31 (157) 
and the location of the minimum of the dispersion curve 



£o 



= 2 + 0(e) « 2. 



(158) 



Although we have only the first two terms in the e ex- 
pansion, these extrapolated values give reasonable re- 
sults consistent with the Monte Carlo simulations or the 
experimental measurements. Furthermore, the next-to- 
leading-order corrections are not too large even at e = 1, 
which suggests that the picture of the unitary Fermi 
gas as a collection of weakly interacting fermionic and 
bosonic quasiparticles may be a useful starting point even 
in three spatial dimensions. 

We have also formulated the systematic expansion for 
the unitary Fermi gas around two spatial dimensions. We 
used the results around d = 2 as boundary conditions 
which should be satisfied by the series summations of 
the expansion over e = 4 — d. The simple Borel-Pade 
approximations connecting the two expansions yielded 
^ = 0.378±0.013 at d = 3, which is small compared to the 
naive extrapolation of the e expansion (Fig. 15). In order 
for the accurate determination of ^ at d = 3, the precise 
knowledge on the large-order behavior of the expansion 
around four spatial dimensions as well as the calcula- 
tion of higher order corrections would be desirable. Once 
this information becomes available, a conformal mapping 
technique, if applicable, will further improve the series 
summations [33]. 



The phase structure of the polarized Fermi gas in the 
unitary regime has been studied based on the e expan- 
sion. We found the splitting point, where the minimum 
of the fermion dispersion curve sits exactly at zero mo- 
mentum, located on the BEG side of the unitarity point 
where akp « 1 and fi/sp ~ —0.5. Then the gapless 
superfluid phase and the superfluid phase with the spa- 
tially varying condensate were shown to exist between the 
gapped superfluid phase and the polarized normal phase 
in a certain range of the binding energy. Our study gives 
a microscopic foundation to the phase structure around 
the splitting point which has been proposed using the 
effective field theory [26]. Moreover, we found the gap- 
less phase with the spatially varying condensate around 
d = 4 exists only at 0.494 < (afcp)~^ < 1 and termi- 
nates near the unitarity limit in contrast to the conjec- 
tured phase diagram where that phase is connected to 
the FFLO phase in the BCS hmit [26]. Our result sug- 
gests that the superfluid phase with the spatially varying 
condensate existing in the unitary regime may be sepa- 
rated from the FFLO phase in the BCS regime. Further 
study will be worthwhile to confirm this possibility. 

We believe that the e expansion is not only theoreti- 
cally interesting but also provides us an useful analytical 
tool to investigate the properties of the unitary Fermi gas. 
As far as we know, this is the only systematic expansion 
for the unitary Fermi gas at zero temperature that exists 
at this moment. The application of the e expansion to 
study the thermodynamics of the unitary Fermi gas at 
finite temperature will be reported elsewhere [34]. 

Note added. Renormalization group argument to 
support our results appeared after we submitted our 
manuscript. Reference [35] studied a renormalization 
group flow equation and found the flxed point describing 
the unitarity limit to approach the noninteracting Gaus- 
sian fixed point in the limit d = 4 or d = 2. Therefore, 
slightly below d = 4 or above d = 2, the finite density 
system associated with the fixed point is weakly coupled. 
The small effective coupling in our Eq. (15) or Eq. (18) 
was correctly reproduced by the renormalization group 
fixed point [35]. 
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APPENDIX: CANCELLATION OF RENORMALONS AT d = 4 - e 

Here we discuss the ultraviolet renormalon [36-38], which is naively present in the diagram depicted in Fig. 14 and 
produces a n! contribution at order e" from the large momentum integration region of a single diagram. The purpose 
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of this Appendix is to show that the n! contribution cancels with subleading contributions from lower-order diagrams 
in the expansion over e = 4 — d. The n + 1-loop diagram in Fig. 14 can be written as 



(A.1) 



where D{p), Ilo{p), and Ila.{p) are, respectively, defined in Eqs. (30), (31), and (37). Introducing these definitions and 
integrating over po, we obtain the following expression for V^+i: 



Vn+1 = 



(-E'fc-f - £fc-f )(^fc+f - gfc+f ) 
4 Jpk Ek_^Ek+E (^Ek-E + Ek+E + f ) 



-E'fe-f + Ek+p + Ji ^Ei_pEij^v 
I {Ei_v + £i_v){Eij^v + ei+p) 



(A.2) 



+ 



{Ei_v - ei_p){Ei_^_P - £j+p) 



ra-l 



Now wc consider the I integration in the bracket. Since the I integration contains a logarithmic divergence at d = 4, 
we subtract and add its divergent piece as 



1 + 



E-, 



1 



•.El 



Ek-R + Ek+i 



+ 1 



E + f / 4Ei_eEi+.{ } 



2ei + Ek_E+E^+E + '-f 



(A.3) 



Ek-E + Ek+B + ^Ji 2ei + Ek-B + Ek+R + f 



The I integration in the second line becomes finite at d = 4 and does not produce singular logarithmic terms. So we 
concentrate on the I integration in the last line, which can be evaluated in the dimensional regularization as 



1 + 
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Ek_z + E^+s. +^ Ji 2£j + Ek-E + Ek+E + -f 
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Then Vn+i in Eq. (A.2) becomes 

^2 r {Ep-ep){E^-e^) 



V, 



n+l 



9_ 

' 4 



EpEq {Ep+Eq+^-^) 



E„ + Ea + ^ 



n n-1 



+ regular terms 



, (A.5) 



where the integration valuables are redefined to be A; — ^ — > p and k + ^ — > q. "Regular terms" coming from the 
second line in Eq. (A.3) vanish in the limit p — !■ cx) or q — > 00. 

If one expanded V"„+i in powers of e and picked up the logarithmic terms, one would find the p integration from 
its large momentum region ^q, £q gives the following contribution to the effective potential at 0(e"): 



4 Jpq EpEq {Ep 



{Ep - e p){Eq - Eg) 

+ Eq+'-^) 



In 



n-l 



6 \ n— 1 



(|)"r(n), 



(A.6) 



which grows as a factorial of n. Whether such a n\ contribution at 0(e") really survives and dominates the large- 
order behavior of the perturbative expansion is a subtle problem. Actually in our case, we can show that the factorial 
contribution of Vn+i at 0(e") cancels because the lower-order diagrams Vi<„ also produce n\ contributions to 0(e") 
with different signs. 

In order to see the cancellation explicitly, we sum up Vn+i in Eq. (A.5) over n to result in 



n=l 



{Ep - ep){Eq - 



E E E 

pq ^p^q 
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+ regular terms 



(A.7) 
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Prom the large momentum region of the p integration, we obtain 

where the n! contribution at 0(e") has been canceled with the subleading contributions from the lower-order diagrams. 

Since the instanton contribution to the large orders, even if it exists, typically has alternating signs, the possibility of 
the Borel summability of the e = 4 — d expansion, which was assumed in Sec. VI, remains open. 
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